Non-negative subtheories and quasiprobability representations of qubits 
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Negativity in a quasiprobability representation is typically interpreted as an indication of non- 
classical behavior. However, this does not preclude states that are non-negative from exhibiting 
phenomena typically associated with quantum mechanics — the single qubit stabilizer states have 
non-negative Wigner functions and yet play a fundamental role in many quantum information tasks. 
We seek to determine what other sets of quantum states and measurements for a qubit can be non- 
negative in a quasiprobability representation, and to identify nontrivial unitary groups that permute 
the states in such a set. These sets of states and measurements are analogous to the single qubit 
stabilizer states. We show that no quasiprobability representation of a qubit can be non-negative for 
more than four bases and that the non-negative bases in any quasiprobability representation must 
satisfy certain symmetry constraints. We provide an exhaustive list of the sets of single qubit bases 
that are non-negative in some quasiprobability representation and are also permuted by a nontrivial 
unitary group. This list includes two families of three bases that both include the single qubit sta- 
bilizer states as a special case and a family of four bases whose symmetry group is the Pauli group. 

.2 

For higher dimensions, we prove that there can be no more than 2 states in non-negative bases 
of a d-dimensional Hilbert space in any quasiprobability representation. Furthermore, these bases 
must satisfy certain symmetry constraints, corresponding to requiring the bases to be sufficiently 
complementary to each other. 

PACS numbers: 03.65.Ta, 03.67.-a, 03.65. Sq 



I. INTRODUCTION 

As an alternative to the standard formulation of quan- 
tum theory in terms of vectors in a Hilbert space, it is 
possible to express quantum states, transformations, and 
measurements as functions on some state space. The 
most common of such representations is the Wigner func- 
tion [T] , which represents the quantum state of a particle 
as a distribution over the classical phase space of the 
particle. However, this function cannot be interpreted 
as a probability distribution as it takes on negative val- 
ues. Such descriptions are referred to as quasiprobability 
representations 2\. 

The occurrence of negative probabilities in the de- 
scription of a quantum state or measurement is often 
thought of as an indication of "quantum-ness" ^ . Con- 
versely, any state or measurement that can be described 
by non-negative (true) probabilities is sometimes said to 
be "classical" . For the Wigner function of particle me- 
chanics, Hudson's theorem shows that a pure quantum 
state has a non-negative Wigner function if and only if it 
is a Gaussian state [H [5] . The discrete Wigner function 
extends these results to finite-dimensional quantum sys- 
tems. For finite odd dimensions, the only pure states with 
non- negative Wigner functions are stabilizer states [51 [7] . 
Negative distributions, on the other hand, are nonclassi- 
cal in the sense that they are contextual [S] and can serve 
as a resource for quantum computation [9^ . 

However, preparing a system in a quantum state de- 
scribed by a non-negative Wigner function is neither a 
necessary or sufficient condition to say that this system 
is classical. It is not necessary, because the Wigner func- 



tion is only one possible quasiprobability representation 
of quantum theory. One can construct a quasiprobability 
representation in which any individual state or measure- 
ment has a non-negative distribution. It is also not suf- 
ficient, as it is possible to make all quantum states have 
non-negative distributions, but this forces the conditional 
probabilities of some measurements to take on negative 
values [12] . To say that a system has a classical de- 
scription, we require the set of preparations, transforma- 
tions, and measurements we are considering to all have 
nonnegative distributions in a quasiprobability represen- 
tation. Consequently, wc aim to determine what subthe- 
ories of quantum mechanics (i.e., theories constructed 
from closed subsets of quantum states, transformations 
and measurements) are non-negative in some quasiprob- 
ability representation, or, alternatively, can be described 
"classically" . 

In this paper, we consider a particular class of subtheo- 
ries containing orthonormal bases both as a set of prepa- 
rations and as a projective measurement. For a fixed 
quasiprobability representation, if the distributions cor- 
responding to the states in an orthonormal basis and the 
projective measurement in the basis are all non-negative, 
we refer to the basis as non-negative. 

The question we are then considering is what sets 
of bases can be simultaneously non-negative in some 
quasiprobability representation. Additionally, we would 
like to determine what groups of unitary transformations 
permute a set of non-negative bases, as it then seems 
likely that circuits composed of such gates can be effi- 
ciently simulated. While our results are primarily ap- 
plicable to single systems, we believe that these results 
may be important in the pursuit of additional classically- 



simulatable subtheories of multiple systems. At present, 
we know of very few nontrivial sets of quantum gates that 
can be efficiently simulated on a classical computer, such 
as the Clifford group and matchgate circuits. Finding 
other subsets of gates that are also efficiently simulat- 
able may provide key insights into the nature of quantum 
computation and the origin of any advantage of quantum 
computation over classical computation. 

While our results are derived in the context of single 
systems, it is important to note that several interesting 
questions remain on the "quantum" nature of single sys- 
tems, as demonstrated by a range of phenomena such as 
the proof of contextuality for a single qubit [TD], mod- 
els of quantum computation that use only one "clean" 
qubit [TT] and recent results on the distillability of magic 
states fSl . We believe that our results show some promise 
in answering such questions. 

This paper is structured as follows. We begin by in- 
troducing quasiprobability representations of quantum 
mechanics in Sec. Ill] and prove some elementary prop- 
erties that any quasiprobability representation of quan- 
tum mechanics must satisfy. We obtain an upper bound 
on the number of non-negative bases of a qubit for any 
quasiprobability representation and relations that any 



non-negative bases of a qubit must satisfy in Sec. Ill A 



We consider examples of ciuasiprobability representations 
that are non-negative for subtheories of a qubit with non- 
trivial transformations in Sec. 



IIIB In Sec. IV 



we gen- 
eralize the theorems in Sec. IIII Al to d dimensions and 
obtain an upper bound on the number of states that are 
elements of non-negative bases for any finite d. We con- 
clude with a discussion on the implications for quantum 
computation in Sec. IV) 



II. ONTOLOGICAL MODELS AND 
QUASIPROBABILITY REPRESENTATIONS 

In this section we motivate the study of quasiprobabil- 
ity representations, and non-negative subtheories within 
them, from a foundational perspective based on ontolog- 
ical models. We present our formalism for quasiprobabil- 
ity representations, and prove some elementary proper- 
ties that any such representation of quantum mechanics 
must satisfy. 

One approach to explaining the predictions of quan- 
tum mechanics is to formulate an ontological model (i.e., 
a hidden variable model) that reproduces some or all of 
the measurement statistics of quantum theory. Such a 
model is defined over some ontic state space A. Prepara- 
tions of quantum states correspond to probability mea- 
sures /x over A. Measurements correspond to sets of con- 
ditional measurement probabilities S,{k\X) of observing 
an outcome k given that a system is in the ontic state 
A € A, respectively. A desirable feature of such a model 



is a revised notion of noncontextuality [^ so that each 
preparation as a density operator corresponds to a sin- 
gle probability distribution and each measurement effect 
corresponds to a unique conditional probability |10| . 

While no noncontextual ontological model can repro- 
duce all of quantum mechanics, it may be possible to de- 
fine a noncontextual ontological model for a subtheory 
Such a subtheory may still capture some of the essen- 
tial phenomena of quantum mechanics. For example, the 
set of states, transformations, and measurements with 
Gaussian Wigner distributions in particle mechanics can 
be described by a noncontextual ontological model [T5] . 
For finite odd dimensions, the single qudit subtheory con- 
sisting of stabilizer states and measurements has an onto- 
logical model, namely the discrete Wigner function [T5] 
restricted to this set [51 [7]. Spekkens' toy theory [T7] 
also provides an ontological model that is in many ways 
analogous to the stabilizer subtheory of a single qubit. 
Despite being "classical", these models allow a variety 
of information processing tasks typically associated with 
quantum mechanics, such as quantum teleportation and 
dense coding; see HH [T7H^ . 

An alternate but related approach to explaining quan- 
tum mechanical predictions is to use a quasiprobabil- 
ity representation. As with an ontological model, a 
quasiprobability representation is defined over a state 
space A (often, but not necessarily, a classical phase 
space) which can be interpreted as an ontic state space, 
and it associates preparations and measurements with 
distributions and conditional distributions, respectively, 
over A. A quasiprobability representation is a faithful 
representation of the density operators and measurement 
effects of quantum theory — that is, the map from op- 
erators on Hilbert space to distributions on A is linear 
and injective. Such representations can maintain non- 
contextuality but nevertheless reproduce the quantum 
predictions because, unlike ontological models, the dis- 
tributions corresponding to such preparations and mea- 
surements are allowed to take on negative values (thus 
the term "quasiprobability"). That is, nonnegativity of 
distributions on ontic states is the classical assumption 
dropped by quasiprobability representations in order to 
reproduce the predictions of quantum theory. 

A quasiprobability representation of quantum mechan- 
ics cannot be non-negative for all preparations and mea- 
surements; this fact is equivalent to the fact that an onto- 
logical model of quantum mechanics cannot be noncon- 
textual [8]. However, a quasiprobability representation 
can be non-negative for a subtheory of quantum me- 
chanics; specifically, the preparations, transformations 
and measurements within the subtheory can all possess 
non-negative probability distributions. In such cases, 
the quasiprobability representation provides a noncon- 



^ While there are ontological models that do not satisfy the revised 
assumption of noncontextuality, e.g., the models in |13lll4j . we 
do not consider them here. 



textual ontological model for this subtheory. That is, 
the existence of a noncontextual ontological model for a 
subtheory and a quasiprobability representation that is 
non- negative for a subtheory are equivalent notions [S]. 
Throughout this paper, wc will restrict our language (for 
the most part) to that of quasiprobability representations 
and the possible existence of non-negative subtheories 
within them, although the reader should keep in mind 
that results for non-negative subtheories in quasiproba- 
bility representations apply identically to a perspective of 
noncontextual ontological models for such subtheories. 

An example of a quasiprobability representation fre- 
quently used in quantum optics is the Wigner function. 
The subtheory of quantum mechanics consisting of Gaus- 
sian states, transformations and measurements, is com- 
pletely described by a noncontextual ontological model 
for which the probability distributions and conditional 
measurement probabilities are all non-negative |15j . This 
subtheory is embedded in a quasiprobability representa- 
tion, describing all possible states, transformations, and 
measurements, but using negative probabilities for non- 
Gaussian ones. 

A quasiprobability representation of quantum mechan- 
ics over a space A is defined by two sets of Hermitian op- 
erators, {f (A)} and {G{X)}, acting on a d-dimensional 
Hilbert space Ud ^,^. The sets {F{X)} and {G(A)} are 
dual frames for the space of operators acting on Hd [5T] ■ 
The quasiprobability distribution associated with a quan- 
tum state p is 



^lpiX) = Tr[pF{\)]eR. 
The support of a state p is the set 

5(p) = {AgA:m,(A)^0}. 



(1) 



(2) 



A point A is compatible with a quantum state p if A G 
S (p) and incompatible with p otherwise. 

For measurements, the conditional quasiprobability of 
an effect E (i.e., an element of a POVM) occurring if the 
system is in the state A is given by an indicator function. 



^E{\)=Tr[EG{X)]e 



(3) 



As a system is always in some state, we require p to be 
normalized, i.e.. 



dXpp{X) — 1 



(4) 



for all states p. Similarly, as some outcome of a measure- 
ment always occurs, we require 



E^^.w = i 



(5) 



for all A e A and all POVMs {Ej}. In order to reproduce 
the Born rule, p and ^ must satisfy 



TripE)^ / dXppiX)CE{X). 

JK 



(6) 



for all states p and ah POVMs {Ej}. 

While negative values of p and ^ are allowed, there 
may be states p or effects E such that Pp{X) > or 
^e{X) € [0, 1] for all A G A, respectively. Such states and 
effects are referred to as non-negative. 

In quantum mechanics, there is a one-to-one correspon- 
dence between orthonormal bases of Hd and projective 
measurements, as any orthonormal basis {p{j) '. j G Z^^} 
corresponds to a set of preparations and to the effects for 
a projective measurement. Motivated by this correspon- 
dence, we define a non-negative basis as an orthonormal 
basis {p{j) : j G Z^} such that for all j G Z^ and A G A, 



CpO)(A)g[0,1] 



(7a) 
(7b) 



That is, each state in a non-negative basis is a non- 
negative state and the projective measurement corre- 
sponding to the non-negative basis is a non-negative mea- 
surement. Considering non- negative bases as correspond- 
ing to both a basis of non-negative states and a non- 
negative projective measurement is the fundamental tool 
that we will exploit to obtain the results of this paper. 
Note that when we consider non-negative bases, we only 
consider orthonormal bases of pure states. 

We are particularly interested in subtheories of quan- 
tum mechanics that contain non-negative bases and a 
nontrivial group of transformations that permute the 
non-negative bases. For qubits, we will establish upper 
bounds on the number of non-negative bases and some 
relations any non-negative bases must satisfy. We will 
also completely classify the possible sets of non-negative 
bases that are closed under a nontrivial unitary group. 
For qudits, we will only undertake the first task (i.e., 
establish upper bounds on the number of non-negative 
bases and relations between any non- negative bases). 

We begin by establishing some properties that any 
quasiprobability representation must satisfy. These prop- 
erties were all proven in Rcf. 10 , but are included 
here for completeness. We first prove that the supports 
S [p{j)] of the states in a non-negative basis {p{j) : j G 
Tjd} must be disjoint. Furthermore, for any A that is 
compatible with one of the basis states, the indicator 
functions £,p(j) (A) are outcome deterministic and corre- 
spond to answering the question "is A compatible with 
P(j)?"- 

Lemma II. 1. Let {p{j) : j G Z^} be a non-negative basis 
of "Hd in a given quasiprobability representation. Then 
the supports {S [p{i)] '■ j G Z^} are disjoint, and for all 
j, k ^ "Ld we have 



epO)(A) = <^,>VAG5[p(fc)] 



(8) 



Proof. Let {p(j) : j G Z^} be a non-negative basis. 
As non-negative bases are orthonormal by definition, 
Tr [p(j)p(fc)] = (5jjt. Consequently, for the preparation 
p{j) followed by a measurement in this basis, the Born 



rule gives 



/. 



S[p(j)] 



d^^J■p{])W^p(k)W ~ 5],k ■ 



(9) 



for all j, k G Z^. As /ip(j)(A) > and is normalized, and 
^p(fc)(A) e [0, 1], the only solutioiFlis as in Eq. ([s]). D 

Quasiprobability representations are convex-linear 
(i.e., they are afSne maps, preserving convex combina- 
tions of states and measurements) |2T]. We now show 
that, for a fixed quasiprobability representation with 
at least one non-negative basis, convex-linearity implies 
that any point A is either compatible with exactly one 
element of each non-negative basis, or is incompatible 
with all elements of every non-negative basis. Any points 
A that are incompatible with all elements of every non- 
negative basis are irrelevant to the description of the non- 
negative subtheory, so we define the space A* C A by 
deleting such points. Furthermore, we also show that 
there exists a unique function g : A* — > M+ such that 
any element that assigns nonzero probability to A as- 
signs probability q{X)- This unique value is given by 
g(A) — d ■ /ii j^(A), that is, d times the probability that 
the maximally- mixed state assigns to the ontic state A. 

Lemma II. 2. For any quasiprobability representation of 
T-Ld in which there is at least one non-negative basis, there 
exists a unique function g : A — >■ ]R+ U {0} such that for 
every non-negative basis {p{j)}, A can be partitioned into 
d -\- 1 disjoint regions {S [p{i)] , Aq : j G Z^} such that 



I^pU)W 



SjmW VA e 5 [p{k)] , 
VA e Ao , 



(10) 



for all i, k G Z^j, where Aq = {A : g(A) = 0}. 



Proof. Let {/o(j)} be an orthonormal basis of T-Ld- Then 
the maximally mixed state can be written as 



-d^-7iT.P^^)- 



(11) 



As the quasiprobability representation is convex-linear, 

'^^iW = ^E/VO)W (12) 

j 

for all A e A. Set q{\) = d ■ //ii(A). 

For any non-negative basis {p(j)}, each term in the 
right-hand side of Eq. ( 12 1 is non-negative by definition 



and so q{\) must be non-negative. Furthermore, any 
point A G A is compatible with at most one element of 
any non-negative basis by Lemma |IL1[ D 



We now show that two quantum states that are ele- 
ments of non-negative bases have disjoint supports if and 
only if they are orthogonal quantum states. 

Lemma II. 3. Let p and lo be elements of (possibly iden- 
tical) non-negative bases in a quasiprobability distribu- 
tion. Then S (p) H S (uj) — ^ if and only if p and uj are 
orthogonal. 

Proof. Let p and w be elements of (possibly identical) 
non-negative bases. By Lemma ILl and IL2 ^p(A) is 



nonzero only for A such that /ip(A) > or q{X) — 0. By 
Lemma n.2 /it^,(A) can only be nonzero for A such that 
q{\) > 0. Therefore, for 



dX p,uj{X)^p{X) = Tr(pw) 



(13) 



to hold, there must exist some A such that Pp{X) ^ and 
A'cj(A) ^ if and only if p and a; are not orthogonal. D 



Lemmas |IL IfjIL 3] provide the basic mathematical tools 
that we use to provide an upper bound on the number 
of non-negative bases of Hd in any quasiprobability rep- 
resentation of quantum mechanics. Any quasiprobabil- 
ity representation over a space A which is non-negative 
when restricted to a subset of bases of Hd must satisfy 
the following: the representation of a non-negative ba- 
sis {p{j)} of Hd must correspond to a partitioning of A 
into d disjoint regions {5[/9(j)]} such that p{j) assigns 
nonzero probability for all A € Aj and probability else- 
where. The measurement in the non-negative basis then 
corresponds to determining which region the ontic state 



III. QUASI-PROBABILITY 
REPRESENTATIONS OF QUBITS 

We now consider the simplest case, namely, quasiprob- 
ability representations of qubits (i.e., two-dimensional 
quantum systems). We will first establish upper bounds 
on the number of non-negative bases for a single qubit 
and necessary relations between sets of non-negative 
bases. We will then construct a complete characteri- 
zation of all the sets of bases that are closed under a 
nontrivial unitary group and are simultaneously non- 
negative in some quasiprobability representation. 

The qubit case allows a geometrical approach, as any 
qubit state p € B{%2) can be written as 



P = 2 (1 + ^(Z') ■ '^) 



(14) 



where r{p) G M'^ is the Bloch vector corresponding to p 
and (T is the vector (A, Y, Z) of Pauli matrices. Orthog- 
onal quantum states then correspond to antipodal Bloch 
vectors, i.e.. 



Here and elsewhere, except possibly on a set of measure zero. 



Tr {puj) = ^ r{p) = -r{uj) 



(15) 



Another useful feature of the qubit case is that any pure 
single qubit state can be uniquely extended to an or- 
thonormal basis of a qubit. Therefore we can represent a 
basis by a single Bloch vector f{j) and denote the basis 
elements by 

{p{j, 7) := ^ [1 + ir{j) • a] : 7 = ±} . (16) 



A. Maximum number of non-negative qubit bases 

We begin by showing that a quasiprobability represen- 
tation cannot be non-negative for more than two bases in 
any plane of the Bloch sphere. That is, for any quasiprob- 
ability distribution, the Bloch vectors corresponding to 
any three distinct non-negative bases (if three such bases 
exist) must be linearly independent. We will then prove 
that if four bases are non-negative in some quasiproba- 
bility representation of a qubit, then they correspond to 
the vertices of a right cuboid. 

Theorem III.l. In any quasiprobability representation 
of a qubit, there are at most two non-negative bases in 
any plane of the Bloch sphere. 

Proof. The proof is by contradiction. Let {p{j, 7) : J = 
1,2,3,7 = i} be three distinct coplanar bases that are 
non-negative in some quasi-probability representation. 
We rotate the Bloch sphere so that the bases are in the 
xz plane of the Bloch sphere and label the bases as shown 
in Fig. [1} 

As the three vectors {f{j) : j = 1,2,3} are linearly 
dependent, there exists a vector r(p) and a,b Q [0, 1] 
such that 

r{p) ^ ar{2) + {1 - a)[-r{2)] 

= br{l) + (1 - 6)f(3) . (17) 

Therefore there exists a mixed state p such that 

p^ap{2,+) + {l-a)p{2,-) 
= 6p(l,+) + (l-&)p(3,+). (18) 

As p is convex-linear, we have 

Pp{X) = a^p(2,+)(A) + (1 - a)Aip(2 -)(A) 

= 6^p(i,+)(A) + (l-5)Aip(3,+)(A). (19) 

for all A e A. 

As the three bases are distinct, neither a or 5 can be 
or 1. As p{l,+) and p{3,+) are not orthogonal, there 
exists A' e 5 [p{l, +)] n S [p(3, +)] by Lemma [lOl Then 
by Lemma II. 2 the first line of Eq. (19 1 gives Pp{X') < 

D 



q{X'), while the second gives Mp(A') = g(A') 



The proof that no three non-negative bases can be 
coplanar in the Bloch sphere rests on the existence of 
the vector in Eq. (17). This vector must always exist 




FIG. 1. Diagrammatic representation of the decomposition 
of the vector r(p) in Eq. (171 in terms of {±f(2)} and 
{f(l),f(3)}. 



vectors in a two-dimensional plane, which would be a 
contradiction. 

Initially, one might expect that this proof could be 
generalized to rule out four non-negative bases in the 
full Bloch sphere (i.e., the existence of four non-negative 
bases would correspond to four linearly independent 
Bloch vectors, contradicting the dimensionality of a 
sphere). This intuition almost always holds. However, 
as we now prove, there is a family of exceptions with a 
high degree of symmetry, namely, if the four bases cor- 
respond to the vertices of a right cuboid (i.e., a solid 
with rectangular faces). Therefore any set of four non- 
negative bases of a qubit must correspond to the vertices 
of a right cuboid, which also implies that there exists 
no quasiprobability representation of a qubit with five or 
more non-negative bases. 

Theorem III. 2. If there are four non-negative bases in a 
quasiprobability representation of a qubit, then the Bloch 
vectors corresponding to these four bases must correspond 
to the vertices of a right cuboid. 

Proof. The proof is by contradiction. Let {p(j, 7) : J = 
1, 2, 3, 7 = ±, 4} be four bases that are non-negative in a 
quasiprobability representation of a qubit. Negating the 
vectors fj as necessary (which corresponds to relabeling 
the basis states), there exists a vector r{p) and a,Sj G 
[0, 1] such that 

r(p) = af(4) + (l-a)[-f(4)] 






(20) 



where X)? *j ~ ^i ^^ otherwise the four vectors {r(j) : 
j = 1, 2, 3, 4} would be linearly independent. 
Therefore there exists a state p such that 



as otherwise there would be three linearly independent 



p = ap(4,+) + (l-a)p(4,-) 
3 
= Ysjp{j,+). 



(21) 



As fj, is convex-linear, 

fip{X) = aAip(4,+)(A) + (1 - a)^p(4 _)(A) 

3 

= E^j'^pO-^+)(^)' (22) 

for all A e A. 

We now use the properties of quasiprobability repre- 
sentations to restrict the values of the coefficients and 



then show that the only solutions to Eq. (201 with the 



appropriate coefficients correspond to the vertices of a 
right cuboid. 

As the four bases are distinct and cannot be coplanar 
by Theorem III.1[ none of the coefficients can be or 1. 



By Lemma II. 2 the first line of Eq. (22) implies fJ.p{X) G 
{ag(A), (1 - a)q{X)} for ah A G A*. 

As p(l, — ) and p{2, — ) are not orthogonal, there exists 
A3 e 5 [p{l, -)] n S [p{2, -)] by Lemma [!L31 Lemma [!L2 
then implies S3 € {a, 1 — a}. Interchanging the roles of 
the first three bases and iterating the above argument 
implies Sj S {a, 1 — a} for j = 1, 2, 3. As ^ ■ Sj — 1 and 
a ^ {0, 1}, we have Si = S2 — s^ — 7^ and a G {\, |} 



Substituting these coefficients back into Eq. (20) and 
swapping the sign of r(4) if a = | gives 



r(4) = f(l) + f(2) + f(3) 

We can rotate the Bloch sphere so that 

ri — (— sin 0,0, cos 0) 
r2 = (sin 0,0, cos 0) 



(23) 



(24) 



for some Be (0, f ). Substituting Eq. (|24|) into Eq. ([23) 
gives ^2,(3) — rj;{4) and fy{3) ~ ry{4). As r(3) and r(4) 
are distinct unit vectors (as they are pure states belong- 
ing to distinct bases by assumption), — ^^(3) = rzi^) = 
cos9. Therefore the vectors {r(l),r(2), — r(3),r(4)} cor- 
respond to the vertices of a rectangle in the plane z = 
cos9, so the vectors {±r(j) : j = 1, 2, 3, 4} correspond to 
the vertices of a right cuboid. D 



Theorems |III.1| and |III.2| show that in any quasiprob- 
ability representation of a qubit, any three non-negative 
bases must have linearly independent Bloch vectors and 
any four non-negative bases must correspond to the ver- 
tices of a right cuboid. We will find that these Theorems 
alone do not completely characterize the sets of bases 
that can simultaneously be non-negative in a quasiprob- 
ability representation. In particular, there are further 
restrictions on the sets of three bases that can simulta- 
neously be non-negative. 



Sets of non-negative bases that allow nontrivial 
transformations 



We now further restrict our consideration to non- 
negative subtheories of a qubit that contain at least 



one nontrivial unitary transformation. By nontrivial, we 
mean that the transformation is not a phase-multiple of 
the identity, which, for a set of two or more bases, implies 
that a nontrivial unitary gives a nontrivial permutation 
of the states in the subtheory. Such non- negative subthe- 
ories are analogous to the subtheory of single qubit sta- 
bilizers states, measurements in the corresponding bases, 
and single qubit Clifford transformations. 

The set of Bloch vectors corresponding to a set of non- 
negative bases can be regarded as a set of pairs of an- 
tipodal points on the surface of a (Bloch) sphere. If 
there are transformations that permute the elements of 
these bases, then they correspond to elements of the point 
group of the set of points on the surface of the sphere, 
which will be a discrete subgroup of 0(3) (except in the 
special case wherein there is only one non-negative ba- 
sis) . For the transformations to have a unitary represen- 
tation, they must have determinant -1-1 [i.e., be elements 
of S0{3)]. As we are interested in subtheories of quan- 
tum mechanics, we restrict to elements of SO{3) rj 

The finite subgroups of SO (3) have been completely 
classified. The only possible groups that can permute a 
set of non- negative bases are [52] : 

• Doo, the group of rotations about an axis and tt- 
fiips around some orthogonal axis; 

• C2 — Z2, the cyclic group of two elements; 

• D2, the symmetry group of a rectangle; 

• D3, the symmetry group of an equilateral triangle; 

• D4, the symmetry group of a square; and 

• Oh, the octahedral group, 

where we have used the results of the previous section 
to ignore groups that are not the point group of a set of 
one, two or three linearly independent vectors and their 
negatives, or the vertices of a right cuboid. 

We now completely characterize the sets of bases that 
are non-negative in some quasiprobability representation 
of a qubit and in addition are closed under a nontrivial 
unitary group. We are explicit with our construction for 
the cases of three and four non-negative bases; quasiprob- 
ability representations that are non-negative for one or 



^ We only consider elements of SO{3) rather than 0(3) because 
we aim to describe quantum mechanics, or a subtheory thereof, 
for which only unitary and not anti-unitary transformations are 
allowed. However, it is possible to extend quantum mechan- 
ics to include anti-unitary transformations, and if one wished 
to consider ontological models describing such an extended the- 
ory, then all transformations in 0(3) (including reflections) can 
be included. Because we are considering subtheories with only 
discrete evolution, the restriction to unitary dynamics is not nec- 
essarily well-motivated. Specifically, because the evolution in a 
discrete theory is not required to be continuously deformable to 
the identity, such an extended model may be reasonable. 




FIG. 2. Bloch vector representation of the transformations Z 
and X that map: a) r(l) <-> -^(2) and — r(l) ^-> — r(2) and b) 
r(l) <->■ — r(2) and — r(l) <-!> r(2), respectively. 



two bases appear as subsets of those for three non- 
negative bases, and therefore we do not consider them 
separately. 



One non-negative basis 

If there is a single non-negative basis (i.e., {|0),|1)} 
up to unitary equivalence), then the group of transfor- 
mations that permute the elements of the non-negative 
basis is the continuous (Lie) group D^ . We can therefore 
find a quasiprobability representation over a space con- 
sisting of two points — a classical bit — for which this basis 
is represented by non-negative probability distributions. 

The transformations of the non-negative basis is gen- 
erated by the group of U{1) rotations about the z-axis 
(which leaves the basis states invariant) and X, which 
corresponds to a bit-flip. 



Two non-negative bases 

If there are two non-negative bases, then the point 
group of the four vertices {±r(l),±r(2)} is generically 
a representation of D2. Expressing the Bloch vectors as 



r(l) = (sin 61,0, cos 6*) , 
r(2) = (-sin6',O,cos0) , 



(25) 



the symmetry group is the group of Pauli matrices. The 
two non-negative bases and the generators of the symme- 
try group are illustrated in Fig. [2] For the special case 
defined by f(l) • r(2) = 0, the four vertices correspond 
to the corners of a square and so the point group is D4, 
which for the above vectors is generated by X and Y^. 



Three non-negative bases 

There are two inequivalent families of sets of three 
bases with nontrivial point groups, D3 and Z2, respec- 
tively. Both families have special cases for which the 



Quantum state 


Support 


P(l,+) 


(+,0), (+,1), (-,2), 


-:3) 


P(l,-) 


(-,0), (-,1), (+,2), 


+,3) 


P(2,+) 


(+,0), (-,1), (+,2), 


-,3) 


P(2,-) 


(-,0), (+,1), (-,2), 


+,3) 


P(3,+) 


(+,0), (-,1), (-,2), 


+,3) 


P(3,-) 


(-,0), (+,1), (+,2), 


-:3) 



TABLE I. List of the six possible supports over the set of 
ontic states {(e, a) : e = ±, a G Z4} for the elements of three 
non-negative bases {p{j, 7) : j = 1, 2, 3, 7 = ±}. 



three bases have a higher degree of symmetry. In partic- 
ular, both families include a unique special case in com- 
mon for which the point group is the octahedral group. 
For this special case, we recover the standard single qubit 
stabilizer states, whose point group is the single qubit 
Clifford group. 

We wish to construct quasiprobability representations 
over some ontic state space A, in which the set of three 
bases {p(i, 7) ■ j — 1,2,3,7 = i} ^re non-negative. By 
Lemma |II.2[ any non-negative basis corresponds to a bi- 
partition of the ontic space A such that each element of 
the basis only has support on one of the partitions. Fur- 
thermore, the probability that any non-negative state p 
assigns to an ontic state A G 5 (p) must be 9(A), which 
is double the probability the maximally mixed state as- 
signs to A. Therefore, with three non-negative bases, A 
can be partitioned into 8 regions corresponding to the 
8 combinations of basis states that are non-negative in 
that region. Without loss of generality, we can integrate 
over these regions and, for convenience, we parametrize 
the ontic state space as the 8 points 



A = {(e, a) : e^±,a e Z4} , 



(26) 



and henceforth we write g(A) — g(e, a). 

Our (arbitrary) choices of six distributions over these 
eight ontic states are listed in Tab. |T] 

The distribution (/(e, a) will have to satisfy certain con- 
straints in order to reproduce the correct quantum me- 
chanical predictions. The distribution must be normal- 
ized over the support of each non- negative state, i.e.. 



51 9(^''^) 

(£,a)e5[p(j,7)] 



(27) 



for 7 = ± and j — 1,2,3. These six equations give only 
four independent constraints. 

In addition, in order for the quasiprobability represen- 
tation to reproduce the quantum probabilities of prepar- 
ing a system in one non-negative basis and then measur- 
ing in another non-negative basis, we also require 

Yl '?(^' a) = 2 [^ + 7i72^(ji) • riJ2)] , 



(e,a)6S[p0i,7i)]n5[p(i2,72)] 



(28) 



for all ji ^ J2 , which are obtained by substituting Eq. ([8]) 
into Eq. ^. Only three of these equations (e.g., ji > J2 
and 7i = 72 = +) give independent constraints. There- 
fore, for the cases we consider, we obtain four indepen- 
dent constraints from Eq. (27 1 and three further inde- 
pendent constraints from Eq. (28). With q{e,a) defined 
on eight points, the seven independent constraints en- 
sure there will be at most a one-parameter family of 
quasiprobability representations for which the desired 
sets of states are all non-negative. 

To find a quasiprobability representation of all states 
and measurements of a qubit that is non-negative for the 
above bases, we need to find operators i^(e, a) and G(e, a) 
satisfying Eqs. ([I]), ^ and ^. As these operators should 
give a non-negative distribution for our chosen bases, we 
require 




f(3) / 



Mp(j,7)(e,a) = Tr[F(e,a)p(j,7)] 

^ U(e,a) if (e, a) e 5 [pO', 7)] 
1 otherwise , 



(29) 



for all e, 7 = ±, J = 1, 2, 3 and a = 1, 2, 3, 4. To find such 
F(e,a) and G{e,a), we find vectors d{a) such that 



d(a) ■ r{j) 



1 if(+,a)e5[p(j,+)] , 
— 1 otherwise. 



(30) 



The operators 



FM = '-^n 



ed{a) ■ a 
G{e, a) = -{1 + ed{a) ■ a 



(31) 



then define a quasiprobability representation that is non- 
negative for the three bases. Note that the operators 
F(e, a) and G(e, a) are identical up to normalization. In 
what follows, we will only explicitly present the operators 
F{e,a). 

We now turn separately to each of the inequivalent 
families and construct quasiprobability representations 
that are non-negative for each of the families. As both 
families include the single qubit stabilizer states as a spe- 
cial case, we will discuss this case after covering the gen- 
eral case for each family. 

Case 1. Up to an overall unitary, the first family 
is the one-parameter family of three bases illustrated in 
Fig. |3] with 



r(l) = (sinfl, 0, cos( 



r(2) 



^(3) = (-isin^ 



sin0, cos6' 



n/3 



,— 2 sin 0, COS ( 



(32) 



for 9 € (0,7r). The corresponding point group is D3, 
which is generated by a ^ rotation about the z-axis 



FIG. 3. lUustration of the Bloch vectors in Eq. (32 1, whose 
point group is D3. 



(which we denote by T) and a tt rotation about the y- 
axis (which we denote by 11). These transformations do 
not commute as D3 is a non-abelian group. 

To define a quasiprobability representation for which 
these three bases are non-negative, we need to find a 



distribution q{€, a) satisfying Eq. (27) and (28). The only 
such distributions are 

q{ + ,0) = qo, 

3 
q{+, o.) ^ - sin^ 6 - 1 + qo , a = 1, 2, 3 , 



9(-,0) 



qo 



■ sin^ . 



g(-, a) = 1 - go - - sin^ 6 , a = 1, 2, 3 , 



(33) 



where qo € [0, 1] is a free parameter. The requirement 
that all probabilities should be in the interval [0, 1] im- 
plies 



< (?o < 2 - - sin^ 



(34) 



which can only be satisfied when sin < |. Therefore 
a quasiprobability representation for which these bases 
are non-negative can only be defined if sin < | . This 



constraint is not a consequence of either Theorem III.l 
or |III.2| and so provides an additional constraint on the 
set of non-negative bases. 

In an ontological model, transformations are funda- 
mentally transformations of ontic states, not epistemic 
states (i.e., states of knowledge). For an ontological 
model of a subtheory of quantum mechanics in which 
(pure) quantum states are epistemic states, this means 
that unitary transformations must supervene on trans- 
formations of ontic states (i.e., must be a consequence 
of some underlying transformation of the ontic states). 
In the models we consider, we will always be able to as- 



a) 

P(l.+) 

c) 
P(2,+) 

e) 

P(3,+) 

g) 
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FIG. 4. (a)-(f) Probability distributions /ip(jy^)(e, a) over 
eight ontic states {(e,a) : e = ±,a € Z4}, for the three bases 
{p{j, "f) ■ j = 1, 2, 3, 7 = ±} with Bloch vectors as in Eq. ( 32 1, 



where the quantum states assign nonzero probabiUty to the 
shaded ontic states. The shadings indicate the two nonzero 
probabilities go (dark grey) and qi (light grey) that are as- 
signed to the different ontic states, (g) and (h) Permutations 
of the ontic states that effect the rotations F and 11, respec- 
tively. 



sume that transformations of ontic states are determin- 
istic (i.e., correspond to a permutation of the ontic state 
space). Conversely, some (but not all) permutations of 
ontic states in these models effect a unitary transforma- 
tion. 

The permutations of ontic states that can effect the 
rotations F and n are illustrated in Fig.l4](g) and (h), re- 
spectively. The permutation in Fig. |4] (hi) only permutes 
the probability distributions for non-negative states if 
^(-l-jO) = g(— ,0), which fixes 



q{e, a) =qo = l- - sin 



g(e,a) 



:= gi, a = 1,2,3, 



(35) 



for e = ±, as illustrated in Fig. H (a)-(f). Note that the 
transformation F always supervenes on the permutation 
in Fig. |4] (g), even if g(+, 0) ^ <z(-, 0). 

To construct a quasiprobability distribution for which 
the three bases are non-negative and described by the 
above probability distributions, we need to find vectors 
d{a) satisfying Eq. ( |30[ ). By examining the permutation 
of ontic states that effects the transformation F, we note 
that d{a) = F^~V(1) for a = 1,2,3 [where Tr is the 
fundamental (spin-1) representation of F], so we need 
only find d(0) and d(l), which can easily be determined 
to be 



d(0) = (0,0, sec 6*) , 

d(l) = (|csc6i,0,-isec6i) 



(36) 
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FIG. 5. Illustration of the Bloch vectors in Eq. (39 1, whose 
point group is Z2. 



The operators F(e, a) are then 



90 



F{e,0)^^il + jseceZ), 



F(e,a) = %F?,-Ml 



47 



7 



CSC ex - - sec ez r}-" 



(37) 



for 7 = ± and a — 1, 2, 3, where 




(38) 



is the spin-i representation of F. The operators in 
Eq. ( 37 1 define a quasiprobability representation that re- 



produces all of quantum mechanics for the qubit and 
is non-negative for the states and measurements cor- 
responding to the three bases with Bloch vectors in 
Eq. (pi. 



Case 2. Up to an overall unitary, the second family 
of three bases have Bloch vectors 



f(l) = (sin 6*, 0, cos 61) , 
f(2) = (-sin6',0,cos6') 
r (3) = (cos (f>, sin (j>, 0) , 



(39) 



for 9 e (0, f ) and cj) € (0, tt), as illustrated in Fig. [5] For 
(/) 7^ "I, the point group of these bases is Z2 (i.e., the 
only nontrivial transformation is a tt rotation about the 
z-axis) . li (f) = ^ , then the point group is the set of Pauli 
matrices, which is a projective representation of D2 (or 
the octahedral group ii 9 = j). 



The only distributions q{e,a) satisfying Eq. (27) and 
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(281 are 



q{+,0)=q{+,3) ■.= qo 
'?(+,2) 

'?(-,2) 



9o - 


— cos 29 cos (h sin 9 

2 2 ^ 


9o - 


— cos 20 H — cos (h sin 9 
2 2 


qi- 


,3)^--qo + -cos29 


1 

9 " 


Qo cos (p sin 9 , 


1 

2 '^ 


" 90 H — cos (j) sin , 



(40) 



where go G [Oi 1] is a free parameter. Requiring all prob- 
abilities to be in the interval [0, 1], we obtain 



<(?(+, l) + g(-,l) =sin6l(sin6'-cos0) < 2, 
< (7(+,2) + g(-,2) =sin6l(sin6' + cos(/;') < 2. 



(41) 



Therefore we require | cos 0| < sin 9 in order for the bases 



in Eq. ( 39 ) to be non-negative in a quasiprobability rep- 
resentation. This is another example of a more restric- 
tive condition on non-negative bases than either Theo- 
rem IIILT] or |IIL2] This condition is also sufficient, as for 



such 0, 0, we can set the distribution to 



90 := 9(7,0) = 9(7,3) 



cos 



91 := 9(7,1) 

92 := 9(7,2) 



smf 



sm( 



sm f — cos ( 



- (sin 9 + cos ( 



(42) 



for 7 = ±. 



As with the models for the bases in Ec^. (32), we can 



view the unitary transformations of quantum states as 
supervening on permutations of ontic states. The distri- 



bution in Eq. ( 42 ) captures the symmetry of the bases 



and the possible transformations, as unitary transforma- 
tions can only supervene on permutations of the ontic 
states for values of 0, (p for which the bases have the ap- 
propriate symmetry. 

For ^ ^, the only nontrivial unitary transforma- 
tion that permutes non-negative states is Z (i.e., a 7r-flip 
about the 2:-axis) . The permutation of ontic states that 
effects Z is depicted in Fig. [6] (g). 

For = ^ , the symmetry group is the group of Pauli 
matrices, generated by X and Z, which is a spin-^ rep- 
resentation of D2. The permutation of ontic states that 
effects X is depicted in Fig. [6] (h). This permutation 
maps (+, 1) to (— , 2), and so we require g( -I- , 1) = q(— ,2), 
which is satisfied if and only if ( 
in Eq. (|42|. 



I in the distribution 



To construct a quasiprobability distribution for which 
these bases are non-negative, we need to find vectors d{a) 



a) 

P(l,+) 

c) 

P(2,+) 

e) 

/^(3,+) 

g) 
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FIG. 6. (a)-(f) Probability distributions ^p(j,y)(e, a) over 
eight ontic states {(e, a) : e = ±, a £ Z4}, for the three bases 



with Bloch vectors as in Eq. (39 1. The shadings indicate the 
three nonzero probabilities, go (checkered), q\ ( light grey) 
and g2 (dark grey), that are assigned to the different ontic 
states, (g) and (h) Permutations of the ontic states that ef- 
fects the unitary transformations Z and X, respectively. Note 
that the permutation in (h) is only valid for the special case 
where the bases have D2 symmetry (i.e., when q\ = 52)- 



satisfying Eq. ( 30 ) . One such set of vectors is 



(43) 



d(0) = (0, CSC 0, sec 6*) , 
d(\) — (esc 6', — CSC (p — cos (p CSC 9, 0) , 
d(2) = (— CSC 9, — CSC (p + cos (/) esc 61, 0) 
d(3) — (O,csc0, — sec6') . 

The corresponding operators -^(7, fc) are then 
i^(7,0) = y (l+7csc0y-|-7sec6'Z) , 
i^(7, 1) = — [1 + 7 CSC 9X — 7(csc (p + cos (p esc 9)Y\ , 

92 

i^(7, 2) = — [1 — 7 CSC 9X — 7(csc (j) — cos (j) esc 9)Y] , 



^(7,3) = y (l+7CSC0y-7sec6'Z) , 



(44) 



for 7 = ±. These operators define a quasiprobability 
representation that reproduces all of quantum mechanics 
for the qubit and is non-negative for the states and mea- 
surements corresponding to the three bases with Bloch 
vectors in Eq. (39). 



Stabilizer states. We now consider the common spe- 
cial case of both families of three bases. For this special 
case, the three bases are equivalent to the vertices of 
a regular octahedron (i.e., single qubit stabilizer states) 
and so have a larger symmetry group, namely, the oc- 
tahedral group (equivalent to the single qubit Clifford 
group). In the appropriate special cases, the distribu- 



tions q{e,a) in Eq. (35) (sin 9 
and 9 



f) andEq. (42) 



) become uniform and g(e, a) = j for all e, a. 
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FIG. 7. (a)-(f) Probability distributions fip{e,a) over eight 
ontic states {{e,a) : e — ±,a £ Z4} for the single qubit sta- 
bilizer states. The shading indicates the ontic states that are 
assigned probability | for a given stabilizer state, (g) and 
(h) Permutations of the ontic states that effect the unitary 
transformations H and P, respectively. 



The single qubit Clifford group is generated by the 
Hadamard gate H and the phase gate P. These trans- 
formations supervene on the permutations of ontic states 
shown in Fig. [t] in the basis r(l) = a;, r(2) = y and 
r(3) = z. Thus, all Clifford transformations can be 
viewed as supervening on permutations of ontic states 
in these models. 



General features of quasiprobabihty representations with 
three non-negative bases 

We now briefly discuss some of the features of these 
quasiprobabihty representations, and how they relate to 
Spekkens' toy theory |T7j and the discrete Wigner func- 
tion [16|. 

Studying non-negative bases that are closed under a 
group of nontrivial unitary transformations has shown 
that Theorems |II1.1| and |1II.2| do not completely char- 
acterize the possible sets of single qubit bases that can 
be non-negative in some quasiprobabihty distribution, 
as neither theorem excludes the bases in Eq. ( 32 1 for 

for 1 > 



16*, for 



sin''6' > I or in Eq. ([3y|) tor 1 > |cos0| > smf 
which the bases are "close" to coplanar. This suggests 
that models for sets of three bases only exist when they 
are sufficiently far from being coplanar. However, as an 



exception to this, the bases in Eq. (32) still admit a non- 
contextual model when sin 6* -^ 0, i.e, when all three 
bases are close to being degenerate. 

Note also that whenever a quasiprobabihty represen- 
tation can be defined such that the three bases are non- 
negative, then any unitary transformations that permute 
non-negative states can be interpreted as supervening on 
permutations of the ontic states. 

For the remainder of this discussion, we focus on our 



quasiprobabihty representation of the stabilizer states. 
Our quasiprobabihty representation of stabilizer states is 
related to the standard definition of the discrete Wigner 
function, though our representation is defined over eight 
points rather than four, as follows. Note that in our 
quasiprobabihty representation, each non-negative state 
is uniquely defined by its distribution over the reduced 
phase space {{+,a) : a G Z4}. Over this reduced phase 
space, the operators in either Eq. ( 37 ) or Eq. ( 44 ) can be 
written as 



F{+,a) = c^p(j,7j,a)-dl. 



(45) 



for a e Z4, where 7j_a is chosen such that (+,a) e 
S [p{j,Jj,a)]- By ignoring the points outside the reduced 
phase space and doing suitable renormalizations (i.e., set- 
ting c = d = i), we recover a typical definition of a dis- 
crete Wigner function [T^] . Note that in this setting, the 
Clifford transformations permute the non-negative basis, 
but do not supervene on permutations of ontic states. 

Our non-negative quasiprobabihty representation of 
preparations and measurements in the stabilizer bases, 
together with Clifford transformations amongst these 
states, also recovers those of Spekkens' toy theory [17] 
when restricted to this reduced phase space. In Spekkens' 
theory, one cannot define suitable ontic transformations 
on which all single qubit Clifford transformations super- 
vene. For our theory restricted to the reduced phase 
space, this is also the case. For example, the permutation 
of ontic states that effects 11 [depicted in Fig. [41(h)] mixes 
the values of e and so cannot be defined on this reduced 
ontological space. Furthermore, anti-unitary transforma- 
tions can supervene on permutations of ontic states, as 
is the case in the toy theory. 

Unlike the toy theory, our model of the stabilizer states 
over the full phase space of 8 points allows all Clifford 
transformations to supervene on permutations of ontic 
states. One can ask the additional question of whether 
all permutations of ontic states that permute the non- 
negative states are allowed in the theory. If one requires 
that only unitary transformations supervene on ontic per- 
mutations, then the answer is no. However, if one al- 
lows antiunitary transformations, provided they permute 
the non-negative states, then all unitary and antiunitary 
transformations that permute the stabilizer states super- 
vene on permutations of the ontic states. In particular, 
the permutation that maps e — > — e and leaves a un- 
changed effects a universal NOT gate. 



Four non-negative bases 

The Bloch vectors 

r(±,±,±) :— (±cos0sin6', ±sin(/!)sin6',±cos6') (46) 

for d,(f) G [0, ^], corresponding to the vertices of a right 
cuboid give a set of four bases. Theorem III. 2 shows 
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that any set of four non-negative bases in an arbitrary 
quasiprobability representation must be of this form (up 
to unitary equivalence). 

The high degree of symmetry makes it easy to con- 
struct a quasiprobabihty representation that is non- 
negative for the above bases for any values oi 9,(j) £ 
(0,^). By examining Eq. (30), it is apparent that up 
to an overall permutation of A, we have 



(47) 



'+++ 



where the Ca are the canonical basis vectors of M'^ . 

We therefore consider a quasiprobability representa- 
tion over six points, {(e, a) : e — ±,a — 1,2,3}. Solving 
Eqs. (27) and (28) for four non-negative bases gives 



9(e,l) 



qi 



1 + sin" 9 cos 2(f> — cos" 9 , 



q{e, 2) =: (72 = 1 - sin^ 9 cos 2^ 
g(e,3) =: 93 = l + cos20. 



cos 



(48) 



for e = ±. Unlike the quasiprobability representation for 
three non-negative bases, this quasiprobability represen- 
tation is valid for all 6*, S (0, §), even when the four 
non-negative bases are arbitrarily close to being copla- 
nar. 

The point group of the four bases for general 9 and 
4> is the group of Pauli matrices. The supports of the 
non-negative states and the permutations of ontic states 
corresponding to X and Z are illustrated in Fig.[8](a)-(j). 

There are three equivalent cases with slightly higher 
symmetry, corresponding to when two parallel faces of 
the cuboid are squares. For example, the face in the z 
plane is a square when </> = ^. In this case, the bases 
are also invariant under the phase gate, P, as shown in 
Fig. 1(1). 

In the special case of the cube (which occurs when 
(/) = J, cos 6* = 4=), the symmetry group of the four 
bases is the octahedral group. The permutations of ontic 
states corresponding to two generators of the single qubit 
Clifford group (the Hadamard H and P) are illustrated 
in Fig. [s] (k) and (1) respectively. 



Note that the bases in Eq. (32) with sm^ 9 = | (i.e., 



cos = ± I ) are three of the four bases corresponding to 
the vertices of the cube, so in this limiting case, a fourth 
non-negative basis (with Bloch vectors parallel to the z- 
axis) can be added. For the model in Eq. (35 1 in the 



limiting case sin'^ 9 — ^, q^ — Q and so only six points 
are assigned nonzero probability. If the points (±,0), 
to which any non-negative state assigns zero probability, 
are ignored, then the distributions in Fig. [4] (a)-(f) are 
identical (up to a relabeling of the states) to those in 
Fig.|8](c)-(h). 



a) 

c) 

p(+,+,-) 



1 2 3 

www 



b) 




fA-, 


-,-) 


d) 




P{-, 


--+) 



2 3 



2 3 



e) 

g) 12 3 

p(-,+,+) 



2 3 

www 
www 
www 



1 2 3 



i) 12 3 

X 



p{-,+,-) 

h) 
P(+,-,-) 

i) 






2 3 



1 2 3 



k) 12 3 

H 



1 2 3 



* i * 


c T , 



^g 



FIG. 8. (a)-(h) Probability distributions over 6 points {(7, k) : 
y — iiz,k — 1, . . . , 3], for the eight quantum states with Bloch 
vectors as in Eq. ( 46 1 corresponding to the vertices of a right 



cuboid. The shadings indicate the three nonzero probabilities, 
gi (light grey), §2 (dark grey) and 53 (checkered), that are 
assigned to the different ontic states, (i)-(l) Permutation of 
ontic states that effect the transformations (i) X, (j) Z, (k) 
H, and (1) P respectively. Note that the permutations in (k) 
and (1) are only valid for special cases with higher symmetry, 
namely, when qi = qs and qi — 52 respectively. 



Example of bases that do not admit a preparation 
noncontextual model: Icosahedron 



We have provided an exhaustive list of sets of bases 
with nontrivial point groups that are non-negative in 
some quasiprobability distribution and noticed that the 
single qubit Clifford group (i.e., the symmetry group of 
a regular octahedron) plays a special role. 

The octahedral group is one of two groups that per- 
mute pairs of antipodal points that are the vertices of 
a Platonic solid, the other being the icosahedral group. 
While Theorems |III.l| and |III.2| show that the pairs of an- 
tipodal vertices corresponding to either the icosahedron 
(or its dual, the dodecahedron) cannot all correspond to 
non-negative bases, we provide an explicit proof of this 
fact, which is illustrative of the general case. 

The idea behind the proof is to use the Bloch vectors 
corresponding to the vertices of an icosahedron to obtain 
two decompositions of a mixed state p in terms of two 
different sets of states and then use Lemmas III.2I - III.3] to 
show that in any quasiprobability representation, some 
of the quasi-probabilities must be negative. 

The twelve vertices of an icosahedron are given by the 
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six vectors 



r(0, a) 



r{l,a) 



V^ + v' 
1 


1 



(!,«</', 0) 



(0,1, a^) 



simple map from 'Hd to an intuitive geometric space like 
the Bloch sphere. Therefore we rely on abstract algebraic 
tools rather than geometric ones. We note that this sec- 
tion is quite technical relative to previous sections of this 
paper, and can be skipped upon first reading. 



^I2,a) = ^==(a(p,0,l) 



for a = ± and their negatives, where (p — {1 + v5)/2 is 
the golden ratio. These vertices satisfy the relation 

(1 - 2a)f(l, +) = af(0, +) - ar(0, -) - (1 - 2a)f(l, -) , 

(50) 

where a = (2 + ^)^^. Therefore we have 

p = ap{Q, +, +) + ap{0, -, -) + (1 - 2a)p(l, -, -) 
= (l-a)p(l,+,+) + ap(l,+,-), (51) 

where p{j,a,p) — ^(1 + prj.a)- For a quasiprobability 
distribution, this implies 

/ip(A) = a/ip(o,+^+)(A) + a^p(o _ _)(A) + (1 - 2a)^p(i _ _)(A) 
= (1 - a)^p(i,+^+)(A) + a^p(i,+,_)(A) , (52) 

for all A G A. Assume that all the bases correspond- 
ing to the vertices of an icosahcdron are non-negative. 
By Lemma |11.3 there e xists a A' G 5[p(l,+,— )] n 
S[p{l,—,—)\. By Lemma [ll.2[ for this value of A', the 
first line is at least (1 — 2a)q{X) > 0, while the second is 
exactly aq{X') > 0. As (1 — 2a) > a, this is a contradic- 
tion. Therefore not all the bases can be non-negative. 



Mg^ A. Linear algebra and convex geometry 

We now introduce some tools from linear algebra. 
These tools enable us to deduce the existence of mixed 
states with multiple decompositions in terms of a set of 
non-negative bases if the bases are related in a particular 
way. Such states play an analogous role to the states in 
Eq. ^ and ^ for qubits. 

Two density matrices, p, a; G B{'Hd), are linearly inde- 
pendent (over M) if 



ap + hid — Q 



0, 



(53) 



and are linearly dependent otherwise. Note that we are 
considering linear independence over the operator space 
B{T-Ld), rather than the space of pure states. We intro- 
duce the term disparate to denote a set A = {p{a^i) : 
a G ZtvjJ G Z(i} of N bases of T-Ld such that any set 
A' obtained from A by removing one element from each 
basis and including the identity is a set of linearly inde- 
pendent operators. That is, for all / G Z^ ( where fj 
labels the element of the jth basis that is not in A'), 



a,j\JTtf, 



Pa.jP{a,j) = -.1 ^ Pa,] =b = y{a,j) . (54) 



IV. QUASI-PROBABILITY 
REPRESENTATIONS OF QUDITS 

For qubits, we have shown that no quasiprobability 
distribution can be non-negative for more than four or- 
thonormal bases. Moreover, up to unitary equivalence, 
the only set of four bases that can be non-negative corre- 
sponds to the vertices of a right cuboid. We now turn our 
attention to the general case and establish restrictions on 
the relation between non-negative bases in an arbitrary 
quasiprobability representation of "H^. We will begin by 
constructing a mixed state p that plays an analogous role 
to the states in Eq. ( 18 1 and (pTj) for qubits. We then 
generalize Theorems lll.l and III. 2 to higher dimensions, 
which establishes restrictions on the relation between any 
set of non-negative bases. We conclude this section by 

,2 

obtaining an upper bound of 2 on the number of states 
that are elements of non-negative bases in an arbitrary 
quasiprobability representation. 

For higher dimensions, the relationship between states 
and orthonormal bases is not as straightforward as for 
qubits, where a state can be uniquely extended to an or- 
thonormal basis. Moreover, for arbitrary d, there is no 



The maximum number of disparate bases of Hd is d+l, as 
each basis contributes d — 1 linearly independent density 
operators and including the identity (or the remaining 
element of any of the bases) gives d^ linearly independent 
density operators. 

We denote the span of a set of density matrices A = 
{p{j)} over M by 



TiA)^{J2^jpij)--CjeR}. 



(55) 



The convex hull of a set of density matrices A is the set 
C{A) = {Y^PjPij) : Pj > 0, )^p, = 1} . (56) 



3 



Denote by dC{A) the surface of C{A) with respect to 
T{A). Then, as we now show, if ^ is a set of disparate 
bases, any expansion of any p G dC{A) as a convex combi- 
nation of the elements of A must assign a zero coefficient 
to at least one element of each of the bases in A. 

Lemma IV. 1. Let A — {p(a,j) : j G Zrf,a G Z^v} he a 
set of N disparate bases. Then for all p G dC{A), there 
exist f G Z^ and paj > such that 



a,j\j^fc 



Pa,jPia,j) 



(57) 
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Proof. To prove this lemma, we prove the contrapositive. 
Let p e C{A) be such that ah the coefficients {pa,j '■ j G 
Z^d in some decomposition 



P = ^Pa,jP{aJ) 



(58) 



a, J 



are nonzero for some value of a. For each a G Z^v, we 
can define fa such that 



PaJc = mmp„j 



(59) 



Then we can rewrite p as 

P= Yl iP'^J-Pc'Jc)pi(^J) + ^YPaJc.- (60) 

Therefore we can set all the Pa.f^ to be equal and strictly 
positive without loss of generality. 

We now prove that p is in the interior of C{A) with 
respect to T{A). Let p' G T{A). Then, by definition, 
there exists /3a. j G M with J^a j Pa-i = ^ such that 



If we set 



P' = YPa..]P{a,j) 

a, 3 



5 — min min j3a j < , 



(61) 



(62) 



then for all p' G T{A), the coefficients e/3a,j + (1 — e)Pij 
in the decomposition 

ep' + (1 - e)p = Y.{ePa,j + (1 - e)Po,j)p(a, j) , (63) 



are non-negative whenever 

Pa,i > or e < 



P^j 



Pi,] - Pa, 



Therefore if we set 



e < mm — - 

1,3 Pij - 



then ep + (1 - e)p' e C{A) for ah p' G r{A). 



(64) 



(65) 



D 



Lemma |l V ■ 1 1 shows that any point p on the surface of 
the convex hull of a set of disparate bases can be written 
as 



a-,3\3=^ia 



(66) 



for some / e Z^ and p^.j > 0. Therefore for all pure 
states G T(-4) [which must be on the boundary of 
T(y^)], the lines defined by 






(67) 



for e G [0, 1] must intersect dCiAd) as gl is in the interior 
of C(yl) and (/> is either in dC(A^ or not in C(^). Therefore 
for all pure states (/> G T{A), there exists e > 0, / G Z^ 
and Pcj > such that 



p:==e0+ — -^1 






(68) 



The state defined in Eq. ( 68 ) can be decomposed as a 



convex combination of different bases in multiple ways, 
which allows us to use it to restrict the relation between 
non-negative bases in an analogous way to the states in 
Eq. ra and (pll) for qubits. 



B. Relation betw^een non- negative bases for qudits 



We now use the state p defined in Eq. ( 68 ) to generalize 



Theorem III.l to qudits. 

Theorem IV. 2. In an arbitrary quasiprobability repre- 
sentation, any set of three mutually non- orthogonal non- 
negative bases must be disparate. 

Proof. Let {p{a,j) : a G Z3, j G Z^} be a set of three 
mutually non-orthogonal bases that are not disparate and 
are non-negative in some quasiprobability representation. 
As the three bases are not disparate, we can relabel the 
elements of {p{2,j) : j G Z<j} such that p(2, 0) G T{A), 
where A = {p{a,j) : a G Z2, j G Z^} is the set of the 
elements of the other two bases. 



Therefore we can use the decomposition in Eq. ( 68 ) 



and relabel the elements of the first two bases such that 
/= (0,0) to obtain 

1 d 

ep{2,0) + —^l^Y.T.Po,jP(^'j) (69) 

a=0 j = l 



for some e G (0, 1). As p is convex-linear, we have 

1 d 



eMp(2.o)(A) ■ 



1 



-qW = Yl X1p"j^p("j)('^) 

a=Oj=l 



(70) 



for all A G A. As p(0, 0) and p(l,0) are not orthogonal, 
there exists A' G 5 [p(0, 0)]n5 |p(l, 0)] by Lemma|!o| At 
this value of A', the right-hand side of Eq. (JTOl) is while 
the left-hand side is at least (1 — e)g(A') > 0, yielding a 
contradiction. D 

Theorem IV. 3. In an arbitrary quasiprobability rep- 
resentation, any set of 4 mutually non-orthogonal non- 
negative bases {p{a,j) : a G Z4, j G Z^} must either: 

• be disparate; or 

• satisfy the relationship 

1-e 



ep(3,0) 



a=0 3 = 1 



P{a,j) 



(71) 



for e — j^ (up to a relabeling of basis states). 
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Proof. Let {p{a,j) : a € Z4,j G Z^} be a set of four 
mutually non-orthogonal bases that are not disparate and 
are non-negative in some quasiprobability representation. 

As the bases are not disparate, then, relabeling the 
bases as necessary, we have p(3,0) € T{A), where A — 
{p{a,j) : a £ 1^3, j G Z^}, that is, the set of the ele- 
ments of the first three bases, which must be disparate 
by Theorem |TV]2] 

Therefore we can use the decomposition in Eq. (68) 
and relabel the bases such that / — (0, 0, 0) to obtain 



ep(3,0) 



2 d 
a=Oj=l 



(72) 



Lemma |II.3 
Eq. ([73| is - 



As /i is convex-linear, we have 

_. 2 d 

Q=0 j=l 

for all A G A. 

As p{3, 1) and p{a,j) are not orthogonal for any a 7^ 3 
and an y j, t here exists A' € S[p{3,l)] C] S[p{a,j)] by 
At this value of A', the left-hand side of 
'^(/(A') > 0, while the right-hand side is at 
least pajq{X') > 0, so paj < ^-^ for all a,j. 

As p(3, 0) and p(a, 0) are not orthogonal for any a 7^ 3, 
there exists A" £ S [p(3, 0)] n S [p{a, 0)] by Lemma IL3 
At this value of A", the left-hand side of Eq. (731 is 
(i^ -I- e)q{X") > 0, while the right-hand side is at most 
i=^2g(A") > 0. Therefore e < ^. 

We now want to show that for each a there must exist 
a ka such that Pa,ka, — "ir- To do this, let {a, /3, 7} be a 
permutation of {0, 1,2}. Then as p{l3, 0) and p{-f, 0) are 
not orthogonal, there exists Xa £ S 
At Aq, the right-hand side of Eq. ([73|) is Pa,k. 
■^-j^q{Xa) for some value of ka. The left-hand side is at 



p(_/3,0)]n5[p(7,0)]. 

q{Xa) 



< 



-7^q{Xa), so we have that there exists ka such that 



least 

Pa,kc = ^^- By considering all permutations of {0, 1,2}, 

we see that this holds for all a. 

We can now show use this same approach to show that 
e = j^. Let {a,/3,7} be a permutation of {0,1,2}. 
Then as p{l3,kp) and p{'y,k^) are not orthogonal, there 
exists X'a £ S [p(/3, fc/3)] n S [^(7, k^)]. At A'„, the right- 



hand side of Eq. ( 73 ) is at least ^^2q{X'a)- The only way 



the right-hand side can be equal to the left-hand side is 
if X'a £ S [p{3, 0)] and e = i-f , that is. 



1 

d+l- 



n 



We now show that in any quasiprobability distribution, 
any number A^ of non-negative bases must still satisfy 
a symmetry constraint. However, as N increases, this 
constraint becomes less restrictive. 

Theorem IV. 4. In an arbitrary quasiprobability repre- 
sentation, any set of N > A mutually non-orthogonal 
non-negative bases {p{a,j) : a £ "^NiJ £ 2rf} rnust ei- 
ther: 

• be disparate; or 



satisfy the relationship 



ep(3,0) 



N-2 d 

Y. Yp^^i 

Q=0 j = l 



p{a,i) 



(74) 



for some e < jy_gT^^ (up to a relabeling of basis 
states). 

Proof. Let {p{a,j) : a £ ZnjJ £ I^d} be a set of N 
mutually non-orthogonal bases that are not disparate and 
are non-negative in some quasiprobability representation. 
Relabeling the bases as necessary, we have p{M, 0) £ 
T{A), where M = N-1 and A = {p{a,j) : a £ ZmJ £ 
Zd}, i.e., the set of the elements of the first M bases. 



Therefore we can use the decomposition in Eq. ( 68 ) 



and relabel the bases such that / is the zero vector to 
obtain 



ep{M,0) 



1 






(75) 



As p is convex-linear, we have 



eMp(M.o)(A) + (1 - e)/x 11(A) ^ ^ ^Poj/Xp(„j)(A) 

oGZm j = l 

(76) 

for all A e A. As p{M, 1) and p{a,j) are not orthogonal 
for all a £ "Lm and j £ I^d, there exists X' £ S [p{M, 1)] n 
S [p{a,j)] by Lemma [ll.3 At this value of A', the left- 
hand side of Eq. (76) is ^^q{X') > while the right- 
hand side is at least pa,kq{X'). Therefore pa,k < 
all a £ l^M and k £ l^d- 

As /9(M, 0) and p{a,Q) are not orthogonal for all 
a £ 1m, there exists A" £ S [p{M, 0)] n S [p{a, 0)] by 
At this value of A", the left-hand side of 
^ + e)q{X") > while the right-hand side 

i)i^<z(A"). a 



l-e 
d 



for 



Lemma II. 3 



Eq. (76) is ( 



is at most (Af 



Upper bound on the number of non-negative 
bases for qudits 



Theorems |IV.2f|TV.4| provide strong constraints on the 
relation between any set of non-negative bases in a 
quasiprobability representation. However, it is unclear 
how to use these theorems to obtain an upper bound 
on the number of non-negative bases in a quasiprobabil- 
ity representation. In order to obtain an upper bound 
(which will not be tight), we change tack and exploit the 
fact that for all A S A there exists an operator F{X) 
acting on "H^ such that 



/i,(A) = Tr(pF(A)) 



(77) 



for all p £ B{Hd)- This will enable us to show that there 
are no more than 2 states that are elements of a non- 
negative basis in any quasiprobability distribution, with- 
out requiring that the bases are mutually non-orthogonal. 



16 



To obtain this bound, we note that any density matrix 
can be written as a hnear combination of the F{X) [14], 
so the -F(A) must be a basis for the space of operators 
acting onHd- 

Theorem IV. 5. For any quasiprobability representation 
ofHd, there are no more than 1'^ states that are elements 
of non-negative bases. 

Proof. Let {IIq, : a € ^d'^} be a trace-orthonormal basis 
of BCHd) and let {F(A^) : /3 G Z^s} be a set of cf lin- 
early independent operators, which must exist as the set 
{F(A)} is a basis for the space of operators acting on T-Ld- 
Then, for all p G B{Hd) and /? £ 2^2 , we can write 



aeZ^2 



(78) 



aeZ,2 



where 2 states p, p' e B{Hd) have the same coefficients 
ga for all a if and only ii p = p' . Therefore we can rewrite 
Eq. IttI) as 



Pp 



(^/3) = ^fa,0ga{p)- 



(79) 



As the {F{Xp) : (3 G ^d^} are linearly independent, 
fa^l3 must be invertible. Therefore for any set of val- 
ues {/x(A^) : /3 e I'd^}, there can be at most one state p 
such that Pp{Xf;) — piXp) for all /3 S 2^2. 

From Lemma |II.2[ any state that is an element of a 
non-negative basis can only assign one of two values to 
any point A £ A, namely, or q{X). Therefore there are 
only 2'^ possible sets of values of p, over {A^ ■ P & ^d^} 
that correspond to elements of non-negative bases. D 

The bound on the number of states that are elements of 
a non-negative basis in Theorem |IV. 5 1 is not tight. For ex- 
ample, not all vectors g correspond to a valid density op- 
erator. In particular, no quantum state can have ga = 
for all a. Furthermore, if an element of a non- negative 
basis assigned nonzero probability to more than cP — d+1 
of the points {A/g}, then as the elements of a non- negative 
basis have disjoint support, at least one of the other ele- 
ments of a non-negative basis would have to assign zero 
probability to all of the points {A^} and so would have 
ga — for all a. Therefore all non- negative states assign 
nonzero probability to between 1 and (P — d + 1 of the 
points {A^}. 

However, even accounting for this does not substan- 
tially decrease the upper bound. Furthermore, for qubits, 
we proved in Theorem |III.2| that no more than 8 states 
can be elements of a non-negative basis. The upper 
bound from Theorem |IV.5| is 16 states, and even exclud- 
ing the combinations of p discussed above only reduces 
the upper bound to 14 states. 



V. DISCUSSION AND CONCLUSION 

We have shown that for any quasiprobability represen- 
tation of a qubit, any three non- negative bases cannot 
be coplanar in the Bloch sphere (i.e., they must be dis- 
parate). Moreover, if there are four non-negative bases, 
then they must correspond to the vertices of a right 
cuboid circumscribed by the Bloch sphere. We provided 
an exhaustive list of all "classical" subtheories of a qubit 
that include states, measurements and nontrivial trans- 
formations. These cases revealed several interesting fea- 
tures. Both families of three bases that are permuted 
by a nontrivial unitary group can only be non-negative 
in a quasiprobability representation when they are suf- 
ficiently "far" from being coplanar. However, there is 



an exception to this behavior, as the bases in Eq. (32) 



are non-negative in some quasiprobability representation 
even when 0^0 (i.e., the three bases are almost degen- 
erate). We have also found that whenever a subtheory of 
qubit states and measurements are non-negative in some 
quasiprobability representation, there exists a quasiprob- 
ability representation in which all unitary transforma- 
tions that permute non-negative states correspond to a 
permutation of the ontic states. 

While we have primarily focused on the qubit case, we 
have also shown that the results for qubits directly gen- 
eralize in that any three mutually non-orthogonal non- 
negative bases in a quasiprobability distribution must be 
disparate and any four or more mutually non-orthogonal 
non-negative bases must either by disparate or satisfy a 
symmetry constraint. In this sense, quantum states and 
measurements with a small amount of complementarity 
can be quite difficult to model in a classical theory. In 
addition, we have obtained an upper bound of 2^ on 
the number of states that are elements of a non-negative 
basis. 

We conclude with some discussion of the implications 
of our results for quantum computation, and some fu- 
ture research directions. While our results have been 
presented in the context of single qudits, they are equally 
applicable to multiple qudit systems. Our upper bound 
on the number of non- negative basis states of a qudit, 
although quite loose, suggests that universal quantum 
computation leads to negativity in any quasiprobability 
distribution. This matches the intuition obtained from 
the specific case of the single qudit discrete Wigner func- 
tion [23J. 

While our higher-dimensional results can be applied 
to quasiprobability representations of multiple qubit sys- 
tems, it is not clear how this approach accords with 
classical simulations of quantum systems. In particu- 
lar, multi-qubit stabilizers can be efficiently simulated 
classically (24j and yet do not correspond to a set of non- 
negative bases in any quasiprobability representation. To 
see this, note that stabilizer states and X and Y measure- 
ments (i.e., in bases corresponding to stabilizer states) 
can lead to violations of a Bell inequality [5S] and so 
cannot admit a locally causal model. 
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A natural way of generalizing a quasiprobability rep- 
resentation for a single qubit to one for multiple qubits 
is to take tensor products of the operators {F(A)} and 
{G(A)} that define the single qubit quasiprobability rep- 
resentation via Eq. (fTl) and (Is]). By construction, such a 
quasiprobability representation will be non-negative for 
all tensor products of the single qubit states with non- 
negative distributions, but may also be non-negative for 
other bases that include entangled states. For the en- 
tangled states to be accessible in a classical subtheory, 
there must be some unitary that permutes non-negative 
bases and maps a non-negative product basis to a non- 
negative basis that contains an entangled state. Such a 
unitary can only be viewed as supervening on a permu- 
tation of ontic states (which could always be done for 
a single qubit) if it leaves the set of tensor products of 



the {-F(A)} invariant under conjugation. Unfortunately, 
it is unclear whether such unitaries exist for any set of 
operators that define our single qubit quasiprobability 
representations, although based on the results of [5^ we 
have some evidence to suggest that they do not. 
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